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Abstract. We present low resolution near-infrared spectroscopy of an unbiased sample of 24 ultraluminous infrared 
galaxies (ULIRGs), selected from samples previously observed spectroscopically in the mid-infrared with the Infrared 
Space Observatory (ISO). Qualitatively, the near-infrared spectra resemble those of starbursts. Only in one ULIRG, 
IRAS 041 14-51 17E, do we find spectroscopic evidence for AGN activity. The spectroscopic classification in the near-infrared 
is in very good agreement with the mid-infrared one. For a subset of our sample for which extinction corrections can be 
derived from Paor and Bry, we find rather high Paa luminosities, in accordance with the powering source of these galaxies 
being star formation. [Fell] emission is strong in ULIRGs and may be linked to starburst and superwind activity. Additionally, 
our sample includes two unusual objects. The first, IRAS F00183-7111, exhibits extreme [Fell] emission and the second, 
IRAS F23578-5307, is according to our knowledge one of the most luminous infrared galaxies in H2 rotation-vibration emis- 
sion. 
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1. Introduction 



Ultraluminous infrared galaxies (see lSanders & Mirabel 19961 
for a review) radiate the bulk of their energy in the far-infrared 
as a continuum emitted by warm dust [T ~ 40 - 60 K). The 
dust is heated by ultraviolet radiation either coming from an 
immense starburst in the central region of these objects or from 
an accretion disk surrounding a black hole. Since their discov- 
ery by IRAS there has been an intense debate about the na- 
ture of ULIRGs an d their possible evolutionary connection to 
quasars. FSanders et al. i 1988a) postulated that the dominant en- 
ergy source in a ULIRG is a dust enshrouded quasar appearing 
in the last phase of the merging of two or more gas-rich disk 
galaxies. Although evidence for both starburst and AGN ac- 
tivity in ULIRGs has been found, the question of what gener- 
ally dominates the luminosity has remained largely unresolved, 
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mainly due to observational difficulties associated with the 
large dust obscuration. Since extinction in the optical is about a 
factor of 10 higher than in the near-infrared (e.g., Mathis 1990l) 
and as the latter provides both starburst signat ures and direct 
spectroscopic indicators of AGN activity (e.g.. 'Marc oni et all 
1994) , near-infrared spectroscopy is an important diagnostic 
for these objects. With the advent of sensitive near-infrared 
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iValdes et al. 2005) began to unveil the nature of these most lu- 



minous galaxies in the local universe. Most near-infrared spec- 
tra of ULIRGs were found to be starburst-Uke, but the starburst 
luminosity traced by the observations appeared typically insuf- 
ficient to account for the full ULIRG luminosities. 

With ISO llKessler et alJll996l) it became possible to study 
even more obscured regions by means of sensitive mid -infrared 
spectroscopy. Fine structure lines in 15 ULIRGs ( Genzel et alJ 
[L998) were observed with the ISO short-wavelength grat- 
ing spectrometer. By using the spectrophotometric mode of 
the photometer on board ISO, low resolution spectroscopy of 
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the mid-infrared aromatic features ('PAHs' ) was obtained for 
a sam ple of 62 ULIRGs jLutz etal.lll998t [iiaopoulou et al. 
1 1999 ). Finall y, 16 ULIRGs of p articularly high luminosity were 
obseryed by iTran et alJ (1200 ih with the low resolution spec- 
troscopic mode of the ISOCAM camera. These studies of 78 
ULIRGs, draw n mainly from the complete IRAS 2 Jy sample 
dStrauss et al.lil99a ,1992.) and its extension, the IRAS L2 Jy 
sample dFisher et alJl995l) . find that 70 - 80 % of the ULIRGs 
are dominated by starburst activity but that there is also a trend 
toward s AGN dominance at the highest luminosities. lLutz et alJ 
( 1 19991) compared ISO classifications (starburst or AGN) with 
classifications from optical spectroscopy and found a good 
agreement if optical LINER spectra were assigned to the star- 
burst group, supporting a scenario in which infrared-selected 
LINERs are due to shocks that are probably related to galactic 
superwinds. 

We carried out follow-up low resolution near-infrared spec- 
troscopy for a subset of the ISO ULIRG sample, accessible 
from La Silla, in order to search for obscured AGNs and to 
compare the ISO indicators with those from the near-infrared 
wavelength regime. We used the Paa line which dominates the 
near-infrared spectra of ULIRGs and the coronal [SiVI] line 
as diagnostic tools to search for broad-line components and 
high energy photons, respectively. We investigate the so-called 
'Bry-deficit' quantifying the starburst contribution t o the tota l 
luminosity and first studied by Goldader et al. ( 1995. ll997albl) . 
and discuss the strong [Fell] emission in ULIRGs and its pos- 
sible implications for the nature of infrared-selected LINERs. 
Finally, we discuss two unusual near-infrared ULIRG spectra. 

The structure of this paper is as follows: In section (2) we 
present the sample. In section (3) we report on the observa- 
tions and data reduction, followed by section (4), discussing 
the results of our near-infrared spectroscopy of ULIRGs. 
Section (5) summarizes our main conclusions. We adopt Ho = 
75kms"VMpc and qo - 0.5. 



2. The Sample 

The unbiased I SO sample, compr ising observat ions made 
with ISO-SWS (Genzel et al. 1998), ISOPHOT-S ("Lutz et al." 
1998 ^ .Ri gopoulouetal. 1999) and ISOCAM-CVF ( Tranet aL 
200 1|). was drawn from the complete IRAS 2 Jy sample 
dStrauss et al."l990'. '1992) and its extension, the IRAS 1.2 Jy 
sample ( Fisheiiet _al, .1995). The selected objects had to fulfill 
the following criteria: L40-120 yum > 10"^ L0, equivalent to 
Lir(8 - 1000 yum) > IO12 Lq\ flux Seo > 1.3 Jy, good ISO 
visibility and (for ISOPHOT-S only) redshift below 0.3. The 
sample of 24 ULIRGs presented here is a subset of targets ac- 
cessible for observations from La Silla during our observing 
runs. We did not bias the sample towards starbursts or AGNs. 
Our ULIRGs cover a redshift range 0.043 -0.327. Paa is shifted 
into the K band for afl objects, except for IRAS F00183-71 1 1 
which shows Pa/? shifted into the H band. The sample and its 
basic properties are listed in Tableffl 



3. Observations and Data Reduction 

Near-infrared long-slit spectra of the 24 ULIRGs were obtained 
with SOFI (Moorwood et al. 1998), the camera/spectrometer 
on the ESO NTT 3.5 m telescope at La Silla/Chile, during two 
runs in November 1998 and August 1999. We chose the red 
grism covering a spectral range from 1.53 - 2.52 yum, includ- 
ing the AGN-indicators Pao- and [SiVI] in the K band as well 
as other diagnostic lines. We used a l"x290" long slit giving 
a spectral resolving power R (= A/AA) ~ 600 at 2 /jm. The 
detector used is a Rockwell HgCdTe 1024x1024 Hawaii array 
with 18.5 micron pixels and the pixel scale along the spatial 
axis is 0.292 arcseconds per pixel. 

In single component systems, the slit was oriented along 
the presumed major axis of the ULIRG, as determined from K- 
band images taken by Rigopoulou et al. ( 1999) and Borne et al.| 
(|2000). In addition, a SOFI K-band image was taken of 
IRAS 19458+0944 which is presented in the appendix of this 
paper. In double/multiple component systems, we oriented the 
slit to get spectra of two components simultaneously. Table ^ 
lists the slit position angles used. Spectra were taken in the 
standard way, nodding the telescope ~ 40"along the slit. The 
total integration time was between 20 to 40 minutes. Sky con- 
ditions during both runs were photometric with a seeing of 
0.45-0.8" and 0.8 - 1 .3", respectively. Data re duction was per- 
formed using IRAF- and ECLIPSE software (lDevillardll997h 
in the standard way. Line parameters (flux, center, width, see 
Tables |5]to|3} were measured by fitting a gaussian profile. We 
estimate flux errors less than 25% 

14 of the 24 ULIRGs contain two or more components (see 
also Table If the separation of the components was suffi- 
cient, we extracted spectra for each component, finally obtain- 
ing 34 spectra. 

4. Results 

The spectra for all 24 ULIRGs observed are shown in Fig. ^ 
The hydrogen recombination line Paa is present in all ULIRG 
spectra except F00183-7111 where Paa is shifted outside the 
K band. In most of the spectra, lines of the Brackett series 
such as Bry at 2.166 jim and BrS at 1.945 //m are observed 
as well. Furthermore, in some ULIRGs we detected the helium 
recombination line Hel (2'S-2'P'') 2.058 ;um. In addition to 
the hydrogen recombination lines, most of the ULIRG-spectra 
show strong rotation-vibration transitions of molecular hydro- 
gen, except F00183-7111 because of its redshift. Transitions 
between the vibrational first excited and ground state, 1-0 S(...), 
dominate, suggesting that shocks or dense photo-dominated 
regions ( PDRs) are the main sources of molecular hydrogen 
emission dSternberg & Dalgarndll989l) . In some cases, lines 



' [Perault ( 1987) calculates the infrared luminosity by using the four 
IRAS-bands: L(8- lOOOyura) = 47rD-[1.8x 10-''*(13.48/i2 +5.I6/25 + 
2.58/60 +/,oo)][W]. 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 

' A few ULIRGs of this sample were observed by IValdes et alj 
1 2005) with SOFI in medium resolution mode in the K band. We note 
for some line measurements discrepancies. Possible reasons for that 
are differences in set-up like position angle and weather conditions. 
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Table 1. Sample characteristics 



Object 




z 


Morpho- 


/ fgQ 


optical 


ISO 


RA(2000) 


DEC(2000) 


PA 








loev 




type 


type 






degree 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


F00183-7111 


12.77 


0.327 


s 


0.11 


LINER' 


AGN 


00 20 34 


-70 55 26 


90.0 


00188-0856 


12.31 


0.129 


D 


0.14 


LINERS 


SB 


00 21 26 


-08 39 27 


0.0 


00199-7426 


12.25 


0.096 


D 


0.08 


?- 


7 


00 22 08 


-74 09 26 


-25.5 


01003-2238 


12.19 


0.118 


S 


0.29 


HII' 


SB 


01 02 49 


-22 51 57 


0.0 


01166-0844 


11.99 


0.118 


D 


0.10 


HII' 


7 


01 19 07 


-08 29 1 1 


-65.0 


01298-0744 


12.28 


0.136 


s 


<0.11 


HiP 


SB 


01 32 21 


-07 29 08 


0.0 


01388-4618 


12.03 


0.090 


s 


0.12 


HII* 


SB 


01 40 55 


-46 02 53 


0.0 


02411+0354 


12.15 


0.144 


T 


0.16 


HII' 


SB 


02 43 46 


+04 06 36 


28.0 


03521+0028 


12.46 


0.152 


D 


0.09 


LINER'' 


SB 


03 54 42 


+00 37 03 


80.6 


04063-3236 


11.97 


0.110 


D 


0.08 




SB 


04 08 18 


-32 28 31 


69.0 


04114-5117 


12.16 


0.124 


D 


0.03 




AGN 


04 12 44 


-51 09 40 


40.0 


04384-4848 


12.32 


0.213 


s 


0.07 




SB 


04 39 50 


-48 43 16 


0.0 


06009-7716 


11.95 


0.117 


D 


0.08 




SB 


05 58 37 


-77 16 39 


103.5 


06035-7102 


12.14 


0.079 


D 


0.11 


HII^ ' 


SB 


06 02 54 


-71 03 09 


245.2 


19458+0944 


12.28 


0.100 


s 


<0.07 




SB 


19 48 15 


+09 52 01 


0.0 


20049-7210 


11.93 


0.126 


S 


0.08 


HII'* 


SB 


20 10 27 


-72 01 11 


-13.0 


20100-4156 


12.55 


0.129 


D 


0.07 


HII^ ' 


SB 


20 13 29 


-41 47 34 


41.1 


20446-6218 


12.04 


0.107 


D 


0.08 


HII* 


7 


20 48 44 


-62 07 25 


74.8 


20551-4250 


11.98 


0.043 


S 


0.15 


HII^ ^ 


SB 


20 58 26 


-42 38 57 


0.0 


23128-5919 


11.96 


0.045 


D 


0.15 


HII^ 


SB 


23 15 46 


-59 03 14 


-3.0 


23230-6926 


12.22 


0.106 


D 


0.08 


LINER- 


7 


23 26 03 


-69 10 17 


22.6 


23253-5415 


12.26 


0.129 


S 


0.09 


HII^ 


SB 


23 28 06 


-53 58 31 


-13.0 


23389-6139 


12.10 


0.093 


D 


0.07 


Hipt 


SB 


23 41 43 


-61 22 50 


-8.0 


F23578-5307 


12.11 


0.125 


S 


<0.05 




7 


00 00 24 


-52 50 30 


90.0 



Col. (1) — Name of Object (IRAS ...). 

Col. (2) — Log of Infrared Luminosity Lir in Lq. 

Col. (3) — Redshift z. 

Col. (4) — S=single; D=double; T=triple. 

Col. (5) — Ratio of IRAS flux at 25 and 60 ^lm. 

Col. (6) — Optical type. 

Col. (7) — ISO type; classification based on the PAH feature jRigopoulou et alll999l:lLutz et alll999l:lTran et all200lh . 
Col. (8) — RA(J2000). 
Col. (9) — DEC(J2000). 

Col. (10) — Position Angle PA of the observed ULIRGs. 

lArmus et a l. ( 1989); "Due et al. ( 1997); ' Veilleuxetal. ( 1999 j) ; lAllen et alJiT99ll) ; ISekiguchi & Wolstencrof3<1993h ; iKewlev et ai]j200ll) . 

t These ULIRGs were reclassified by Lutz et alj [l999l) . 

?: Quality of these spectra was not good enough for a classification. 



of [Fell] with rest wavelengths 1.257 and 1.644 //m are ob- 
served in the H and/or K band. In two low redshift ULIRGs 
— IRAS 20551-4250 and IRAS 23128-5919N/S — CO ab- 
sorption bands at 2.29 and 2.32 jum are still present in the part 
of the K band not yet strongly affected by the thermal back- 
ground radiation. These CO absorption bands are produced in 
the atmosphere of red giants and supergiants. 

4. 1 . Search for buried AGNs 

The prime scientific motivation f or our research is the search 
for 'hidden AGN' postulated by ISanders et all (E^SSa). Our 
aim is to compare near-infrared evidence for AGN activ- 
ity with evidence from ISO mid-infrared observations of ou r 
sample ULIRGs (iRigopoulou et all 1 19991: iTran et alJ 1200 ll) . 
Additionally, we can quantify the extent to which near-infrared 
studies are affected by extinction. The near-infrared regime 



contains two indicators linked to the presence of a 'hidden 
AGN' : rSiVI] and the hydr ogen recombination lines (e.g., Paa, 
Pa/I). iMarconi et alJ (Il99 4) showed that a dust enshrouded 
AGN can be detected using the coronal [SiVI] line at 1 .962 fim 
which is observed in most (though not all) local Seyferts. The 
lower ionization potential of this species is 167 eV, therefore, 
only the hard UV-spectrum of an AGN and not a starburst 
can be responsible for th e ionization up to Si^ ^. Spatially re- 
solved observations (e.g.. lMaiolino et al]ll99 8^ and photoion- 
ization models of e.g. ISO-observations jAlexander et al. 1991 
I2OOOI) are consistent with the formation of this line in the in- 
ner part of the narrow-line region. Near-infrared hydrogen re- 
combination lines, especially Paff and Pag, allow a deep search 
for AGN broad-line regions (e.g.. lVeilleuxeTlLl ll997b'l This 
method is obviously sensitive mainly to type 1 AGN while, 
the BLR in type 2s remains usually undetectable out to 4 jum 
jLutz et alJl2002l) . In contrast, [SiVI] emission from the NLR 
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Table 2. Emission-line fluxes 



Object 


[Fell] 


Pa/3 


[Fell] 


Ht 


Paff 


H2 


BiS 


H2 + [SiVI] 


Hi 


Hel 


H2 


Bry 


H2 




1.257 


1.282 


1.644 


1.835 


1.876 


1.891 


1.945 


1.957 


2.033 


2.058 


2.122 


2.166 


2.223 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 



F00I83-7I11 


0.228: 


0.115: 


0.376: 






















00I88-0856N 








0.079: 


0.71 


0.039:: 




0.123 






0.104:: 






00I88-0856S 




























00I99-7426N 










2.07 


0.218 


0.066:: 


0.251: 


0.115:: 




0.462: 






00I99-7426S 










0.126: 


















01003-2238 










0.991 






0.190: 












01I66-0844S 




























01I66-0844W 










0.075 


















01298-0744 










0.072 


0.006: 




0.017 






0.025: 






01388-4618 








0.12: 


1.495 


0.083:: 




0.220 




0.040:: 


0.203 


0.182:: 




0241l+0354Main 










2.906 




0.188 


0.099 




0.213: 


0.125: 


0.408: 




02411+0354W 










0.680 




0.063: 


0.039:: 












03521+0028 








0.162 


0.764 






0.170: 












04063-3236NE 








0.087: 


1.08 




0.055:: 


0.134 






0.123: 


0.104: 




04063-3236SW 










0.207 


















041 14-5 117E 










0.398 






0.089: 


0.053:: 




0.057:: 






041 14-5 117SW 










0.082 


















04384-4848 










1.01 


















06009-7716E 










0.128 


















06009-77 16W 










1.206 




0.102 


0.052 








0.170:: 




06035-7 102W 






0.232: 




1.122 


0.078:: 




0.428 


0.227 




0.480 






06035-7 I02E 










0.231 






0.516: 












19458+0944 






2.079: 


0.911:: 


8.855 




0.460: 


1.276 


0.631: 


0.484: 


1.276: 


1.166: 




20049-7210 










0.172 


















20100-4156* 










0.445 


















20446-6I28N 








0.050 


1.616 


0.040: 


0.099 


0.117 


0.025 


0.899 


0.092: 


0.206 




20446-6I28S 




























20551-4250 






0.257 




1.606! 


0.011! 


0.128! 


0.428! 




0.187 


0.465 


0.199 


0.149 


23I28-5919N 






0.960 




13.06! 




0.695! 


0.441! 




0.802 


0.318 


1.191 


0.417 


23128-5919S 






2.038 




26.33! 




1.578! 


2.696! 


0.451 


1.547 


1.306 


2.675 


0.651: 


23230-6926 










0.221 


1.561 




0.138 


0.278 


0.188 


0.092 


0.247 


0.190: 


23253-5415 








0.075 


0.269 






0.101 






0.175 






23389-6139 






0.600 


0.708: 


2.232 


0.091 


0.121: 


0.653 


0.187: 




0.581 


0.212 


0.341:: 


F23578-5307 








1.395 


2.768 


0.423:: 




1.65 


1.136 


0.966: 


3.735 







Col. (1) — Name of Object. 

Col. (2)-Col. (14) — Line flux in 10-"eig s"' cm--. 

! — Liite in region of poor atmospheric transmis.sion between H and K band. 



can be detected for both type 1 and type 2 AGN. Both meth- 
ods provide a further che ck o f the mid-infrared method used by 
iRigQPOulou et alJ (1 19991) and lTran et alJ (1200 lb which is based 
on the ratio of large scale PAH emission and AGN-heated dust. 

Due to the SOFI spectral resolution of ~ 600, 
[SiVI] is partially blended with the molecular hydro- 
gen line H2V=1-0S(3) at 1.958 fim. Models of flu- 
orescent excitation of molecular hydrogen predict an 
H?, v=l-0S(3)/H7. v=l-0S(l) flux ratio less than unity 
(typically ~ 0.7; iBlack & van Dishoeckl Il987h . while col- 
lisional excitation models predict a ratio that ranges from 



0.5 to 1.4 ( Shull & HoUenbach 1978: Black & van Dishoeck 
ll987HSternberg & Dalgarnolll989t) . Similar to IVeilleux et al 



lll997bl I 1999a) . we have assumed conservatively that 
[SiVI] is present in galaxies where the apparent flux 
ratio of Hz v= 1 -0 S(3)/H2 v= 1 -0 S(l) is larger than 
1.5. IRAS 04114-5117 is a double nucleus system 
which most likely hosts an AGN in one of its two nu- 
clei. In the East component (IRAS 041 14-51 17E), the 
H2 v=l-0S(3)/H2 v=l-OS(l) ratio is 1.6, consistent with 
the presence of an AGN. Also, the observed half width of 
the blend of [SiVI] and H2 v=l-0S(3) lines in that source 
is higher than the width of the H2V=1-0S(1) line alone 
(Table 1^. By fitting multiple Gaussians we have estimated 
the fluxes for the various lines: values of 0.061 x 10"''*erg 
s~' cm"^ and 0.0165 x lO^'^^erg s"' cm"^ were found for 
the molecular hydrogen and [SiVI] lines, respectively. The 
[SiVI] Une flux renders an observed luminosity of 1.4 x 10^ 
L0 for the [SiVI] line. The presence of an AGN in this galaxy 
is also confirmed from the mid-infrared classification using 



the RA^H-diagnostic feigopoulou et alJ Il999l) . The fainter 
component of this system, IRAS 041 14-5 117W, shows no 
evidence of AGN activity. 

For all sources, the hydrogen recombination lines Paa and 
Pa/3, are narrow with an intrinsic full width at half maxi- 
mum (FWHM) of 700 kms^' at most (see Table for flie 
observed FWHM). None of the targets show evidence for 
an additional broad line component (Au a few 1000 km/s). 
In most of our ULIRGs, the Pag line shows a blue wing 
(see also [V eifleux et al."l997bl Il999bl iBurston et al]l200ll) . 
iMurphv et aL U999. .200 L) resolved similar blue wings and 
proposed an identification with two helium recombination lines 
at 1.8686 and 1.8697 u m. Transitions of mole cular hvdrogen 
may contribute as well (iVeilleux et alJ[l9'99bll . Hence, we do 
not consider such asymmetric blue wings as evidence for an 
AGN. 

To better quantify the non-detection of a broad-line com- 
ponent in our sample, we derived 3cr upper limits for the flux 
of the broad-line component of Paa, assuming a line width 
of 3000 kms '. The line width of 3000 kms ' is a typical 
value for Seyfert 1 galaxies and consist ent with the detections 
of broa d Paa lines in some ULIRGs bv IVeilleux etalJ ( Il997bl 
1999b). From this, we determined an upper limit for the ratio 
of the broad line flux and the measured narrow-line flux. A ma- 
jor part of our sample has a limit on this flux ratio lower than 
one (21/29) and almost all (28/29) lower than two. The non- 
detection of broad-line regions is significant considering that 
ratios of broad and narrow flux significantlv above 1 are com- 
mon in Syls (e.g., Stirpe 1990, where the median ratio is 20 
for objects not exceeding a luminosity of M,, - -22 mag), and 
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Table 3. Observed full width at half maximum of emission lines 



5 



Object 


[Fell] 


Pa/3 


[Fell] 


Hi 


Paa 


Hj 


Bi<5 


H2 + [SiVI] 


Hi 


Hel 


H2 


Bry 


H2 




1.257 


1.282 


1.644 


1.835 


1.876 


1.891 


1.945 


1.957 


2.033 


2.058 


2.122 


2.166 


2.223 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 



F00183-7111 


1000:: 


600:: 


1100:: 






















00188-0856N 








572.6: 


481.2 


526.7:: 




509.4 












00188-0856S 




























00199-7426N 










585.3 


717.9 




599.6: 


417.3:: 




823.2: 






00199-7426S 










513.7: 


















01003-2238 










542.1 






667.6: 












01166-0844S 




























01166-0844W 










575.0 


















01298-0744 










634.4 






527.8 






532.5: 






01388-4618 








682.6: 


540.6 


656:: 




607.2 






425.0 






02411+0354Main 










519.3 




494.9 


846.9 




506.3: 


401.3:: 


385: 




02411+0354W 










468.2 




461.8: 


968.0:: 












03521+0028 








433.7 


663.9 






577.5: 












04063-3236NE 








769.6 


587.0 






500.2 






377.9:: 






04063-3236SW 










630 


















041 14-5 117E 










569.8 






945.9: 






634.7 






041 14-5 117SW 










853.7 


















04384-4848 










630.7 


















06009-7716E 










573.8 


















06009-77 16W 










558.2 




598.5 


598.5 








381.0:: 




06035-7 102W 






600.1: 




706.6 


565.9:: 




684.2: 


866.3 




794.1 






06035-7 102E 










592.9 






380.8:: 












19458+0944 






982.9: 


542.4:: 


524.5 




458.6: 


613.6 


730.1: 


407.2: 


561.4: 


578.7: 




20049-7210 










585.2 


















20100-4156 










892.0 


















20446-6128N 








484.4 


502.1 


581.7: 


461.8 


480.0 


551.7 


557.5 


508.8: 


643.4 




20446-6128S 




























20551-4250 






541.0 




582.2! 


526.2! 


481.2! 


592.6! 




523.8 


544.9 


515.9 


453.3 


23128-5919N 






629.8 




598.4! 




568! 


578.4! 




466.0 


336.3 


451.1 


318.7 


23128-5919S 






755.6 




760.3! 




719.4! 


1832! 


614.9 


578.7 


613.6 


655.5 


750.3: 


23230-6926 








863.6 


585.2 




639.9 


714.6 


1152.8 


409.4: 


578.2 


606.3: 




23253-5415 








631.7 


635.7 






609.9 






459.6 






23389-6139 






803.9 


1025.1 


667.9 


350.2 


653.3: 


817.0 


486.2: 




674.5 


336.4 


672.2:: 


F23578-5307 








8824 


623.7 






593.2 


549.9 


706.4: 


713.9 







Col. (I) — Name of Object. 

Col. (2) - Col. (14) — Observed full width at half maximum, in km/s. The instrumental line width is about 500 km/s. 
! — Line in region of poor atmospheric transmission between H and K band. 



have also been observed in a few ULIRGs (IRA S 13451 = 2.3 
IRAS 23499 = 2.06 and IRAS 13305 = 0.86. IVeilleiix etalJ 
Ei223[i999^. IRAS F00183-7111, which was classified as 
AGN by ISO, cannot be classified unambiguously in the near- 
infrared due to its redshift (precluding observation of [SiVI]) 
and faintness (precluding a good measurement or strong limit 
for a broad component to Pa/3). 



Overall, we find evidence for AGN activity in one out 
of the 24 ULIRGs of the present sample. Additionally, there 
is good agreement between near-infrared and mid-infrared 
classifications. None of the 17 ULIRGs classified by ISO as 
starburst show any AGN activity in the near-infrared. Small 
number statistics provided us with only 2 ISO AGN in our 
sample — one (IRAS 041 14-51 17E) also has near-infrared 
evidence for AGN activity, the other (IRAS F00183-7111) 
is uncertain. This is consistent with the good agreement 
found b etween optical an d near-infrared jveilleux et al. 1997 b. 
Il999bl) ISO and optical jLutz et alJll999l) and ISO and hard 
X-ray (" Genzel & Cesar skv' '2000') classifications, suggesting 
that strong AGN in ULIRGs are often identifiable over a 
wide wavelength range. Near-infrared and ISO diagnostics 
also agree in putting most optical LINERs into the star- 
burst category. This is true for three of the four optical 
LINERs observed with SOFI, the fourth is again the am - 
biguous IRAS F00183-7111. As argued by Lutz et al. (1999) 
and Taniguchi et al > 1999.) . this suggests that optical LINERs 
from infrared selected samples are linked to the starburst phe- 
nomenon. We will revisit this issue below when discussing 
[Fell] emission. 



4.2. Starburst luminosity and the 'Bry-Deficit' 

Qualitatively, our results confirm previous findings that the ma- 
jority of near-infrared spectra of ULIRGs are starburst-like. 
iGoldader etal] (Il99.5[ Il997albl) shidied quantitatively in the 
near-infrared whether the starburst activity traced is sufficient 
to power the bolometric luminosity of infrared galaxies. To 
quantify the relative strength of the ionization source, they 
used the Bry-luminosity (which is proportional to the ioniz- 
ing photon rate) in conjunction with the infrared luminosity 
LiR which is assumed to be equal to the bolometric luminos- 
ity of the (U)LIRGs. Bry has been extinction corrected using 
either near-infrared continuum colors or E(B-V) derived from 
the optical hydrogen recombination lines Ha and H/3. Goldader 
et al., found that in luminous infrared galaxies (LIRGs), the 
starburst activity seen in the near-infrared can account for the 
total bolometric luminosity. In ULIRGs, however, they found 
that the starburst activity detected is insufficient to account for 
the bolometric luminosity if similar starburst properties (initial 
mass function, star formation history) to those for the LIRGs 
are assumed. This 'deficit' in Bry luminosity or ionizing lumi- 
nosity, is usually a factor of about 3, although it can reach up 
to 10 in individual cases. A Ukely interpretation is that a signif- 
icant part of the ULIRG's power source has not yet been seen 
in these observations bec ause the central regions are opt ically 
thick in the near-infrared. IVeilleux et all i ll 997hl [l999H) have 
investigated the deficit of the ionizing luminosity using the 
Paff line. They determined extinction (in a foreground screen 
model) using optical/optical and optical/near-infrared hydro- 
gen recombination-line pairs. They confirm the deficit found 
by <Goldader et aL a 99.5. .1997a.bJ for ULIRGs that are cool 
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IRA3FD0ia3-7111 IRASOO iee-0856N 




1,3 1.6 
IRAS00188-0e56S 



1.8 1,4 



1,5 1,6 1.8 1.9 3 8,1 S.S 

IRASOO 199-7436N 




1,6 1,8 1,9 2 2.1 Z.S 1,4 

1RAS00199-7436S 



.1 I.. 

1.8 " i.a 1.9 z 3.1 z.a 1.4 

IRAS01 16B -QB44S 

- o.oa 




1.5 1.6 1.9 

IEAS0116S-0844W 



z. 1 z.a 




1.8 l.i 
IRAS0129B-0744 




' H 0.O4 H 

HO.Ol - 
1.5 1.6 ** 1.! 

IRAS013BB-461B 



Z. 1 2.3 



0,01 yi-'T/v^^'v-wAA-'WVi/'''"''^^^ 







1,6 1,« 1,8 1,9 8 

IRAS024 11 +0354Main 





1.6 IB 1.9 S 3.1 S,S 3.3 

IRA30241H-0354W 
0,2 

0.1 




1,5 1,8 1,9 2 2.1 Z.S 
IRAS035S 1 +002B 



1,6 1,8 1 

IRAS04063-3236NE 




i1l 



1.5 l.B 1.9 Z S.l 

IRAS04063-3236SW 

- I I I 

o.i - 



1.5 1.6 1.9 2 

IRAS041 14-51 17SW 



2.1 2.2 1.4 





1,8 1,9 



1^ 



1,5 1,6 1,7 1,8 1,9 2 



Fig. 1. ULIRG near-infrared spectra - flux density as function of the rest frame wavelength A. For each object we show in 
two panels the data taken in the H and K band and shifted to the rest wavelength. Flux density /a is shown in units of 10 "erg 
s"'cm"^/im ' and rest wavelength A in fim. 



in the IRAS 25 yum/60 jum color but do not find a deficit for 
warm'* (mostly AGN-like) ULIRGs, a category excluded by 
[Goldader et^al. ( 1995, 1997a b) in their initial analysis. Recent 
work by IValdes et alJ ( l2005l) report this deficit using the Paa 
and Bry line as well. 

We used the dominant Paa line to revisit the deficit of the 
ionizing luminosity for our sample of mostly starburst-like and 

" (fv(25 ^m)/f,(60 um')>0.2: ISMiders et alligSSbh 



'cool' ULIRGs. To determine the intrinsic flux of this line, we 
corrected for extinction using the foreground absorber model. 
We adopted the extinction coeffcients of Veilleux et al. ( 1997a) 
and three different line pairs: Ha/HjS, Paa/Ha and Paa/Bry. 
For the calculations we assumed case B recombination, a tem- 

J ierature T ~ 10 "* K and an electron density ~ 10"* cm"'' 
Hummer & Store vl ll987l) . Difi'erent line pairs often provide 
significantly different extinction values. In principle, the ex- 
tinction derived from the purely near-infrared line pair Paa/Bry 
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IRAS0600e-77ieE 
30.06 



IRAS06009-77ieW 

0.3 
0,3 




IS l.« 1.9 S 3.1 

IRAS23 128-56 ISN 

3 



1.5 1.6 



1.7 1.9 8 3.1 s.s a.a S.4 

IRAS33 130-59193 





1,6 l.T 1,8 1,9 S 2,1 3.3 8.3 £.4 

IRAS33353-5415 



1.5 l.B l.B 1.9 2 Z.l 2.2 

IRA3233a9-6139 



1.5 l.B l.B 1.3 2 2.1 Z.a 

IRASP23578-5307 




1.6 l.B 2 2.1 2.3 1.4 



1.6 l.B 2 2.1 2.3 



Fig. 1. (continued). 

is significantly higher than that from the purely optical Ha/Hfi 
or the optical/near-infrared Paa/HQ-, indicating that the near- 
infrared recombination-line emission originates in deeper more 
obscured regions and does not follow a simple screen model. 
Technical reasons like different aperture sizes may also af- 
fect the line ratios derived from independent optical and near- 
infrared observations. On the other hand, Paa and Bry are sep- 
arated by less than 0,3 fim, and thus give useful extinction con- 
straints only for high Ay (>10) because of the small differential 
extinction. For some of our sources, Bry is redshifted into the 
less favourable long end of the K band. Since we are prob- 
ing for obscuration effects that are significant even in the near- 
infrared, and since both hnes are taken from the same spec- 



trum, thus minimizing systematic errors, the Paa/Bry ratio is 
valuable for studying the ionizing photon deficit in ULIRGs, 

In Table0]we summarize the results for those sample galax- 
ies for which we estimated extinctions, and give references to 
the optical spectra used. In Figure |2l we plot the Paor luminosi- 
ties (corrected for extinction based on estimates from the vari- 
ous line pairs) as a function of infrared luminosity. 

Our main results are: For most of the ULIRGs with extinc- 
tion correction based on purely optical Ha/H/3 and optical/near- 
infrared Paa/Ha, we detect a deficit of Paa flux with respect to 
what one would expect f or a starburst (solid Une in Fig.lp , thus 
confirming the results of lGoldader et alJ ( Il995lll997albt) . 
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Table 4. Colour excess E(B-V) as derived from different re- 
combination line pairs 



Object 


E(B-V) 


E(B-V) 


E(B-V) 


Ref. 




Ha 


Pao' 


Paa 






tip 








(I) 


(2) 


(3) 


(4) 






1.46 


1.59 




1 


01003-2238 


0.71 


1.43 




1 


01 166-0844 


0.32 






1 


01298-0744 


0.71 


0.77 




1 


01388-4618 






4.64 




0241 1+0354Main 


0.66 


1.41 


6.37 


1 


0241 1-I-0354W 


0.79 


1.07 




1 


03^21+0028 


2.36 


2.43 




1 


04063-3236NE 






1.85 




04063-3236SW 










060nQ-77 1 6F 










06009-77 16W 






6.40 




0603'S-7102W 


1.36 


0.85 




2 


0603'S-7102F 


0.94 


0.08 




2 


1 945 8 -F 0944 






5.54 




20100-4156 


0.96 


1.75 






20446-6 128N 






5.23 




20446-6 128S 










20551-4250 


1.06 


0.48 


4.88 


2 


23128-5919N 


1.03 


1.45 


1.19 


2 


23128-5919S 


0.76 


1.77 


2.46 


2 


23230-6926 


1.30 


1.28 


4.71 


2 


23389-6139 


1.01 


1.19 


1.67 


1/2 



Col. (5) — Referen ce s for the optical line measurements. 1: 
IVeilleux etaT]h999j) . 2: IDuc et al. 1 1997). 

The calculations are based on case B recombination, a temperature 
T ~ W* K, an electron density «<, ~ lO** cm"^^ I H ummer & Stored 
Il987h and the foreground absorber m odell. We adopted the extinction 
coefficients of lVeilleux et alJ 



For the subsample where the extinction correction came 
from purely near-infrared Paa/Bry line ratio, there is defi- 
nitely no Paa-deficit but rather on average an excess of Paa 
flux. It is difficult to estimate to which extent this excess 
would be reduced if a similar near-infrared extinction correc- 
tion were to be applied to the lower luminosity reference sam- 
ple — the reference objects are typically dusty starbursts for 
which such 'wavelength-dependent' extinction may be present 
as well, at a lower level than in the ULIRGs (see the case 
of M 82 (McLeod et al."l993') which also shows that this 
'wavelength-dependent' extinction is an artifact of enforcing 
a screen model). A tendency may be present in our data for 
galaxies with high Paff/IR to have detections of the Hel line 
at 2.058 jjm (Fig. Despite vario us other comp lex effects 
in the interpretation of this line (e.g.. lShieldslll993l) . this may 
be consistent with younger bursts having both harder ioniz- 
ing radiation and a higher ratio of ionizing to bolometric lu- 
minosity. In a^reern ent with the results from ISO spectroscopy 
(iGenzel et al 1998'), our findings suggest that the deficit in ion- 
izing luminosity derived from near-infrared spectra of ULIRGs 
is indeed an obscuration effect, such that much of the 'missing' 
luminosity is in the form of highly obscured (Av ~ 10-20 
mag) starburst activity. 




11.6 11.8 12.0 12.2 12.-1 12.6 
log Lir/L© 



12.8 



Fig. 2. Logarithm of the Paor-luminosity plotted as a function 
of the logarithm of the far-infrared luminosity Lir. The solid 
line repr esents the relation derived by Goldader et aL ( 1995.. 
Il997albl) for starburst LIRGs, transformed to Paa assuming 
case B ratios. The different symbols represent the different ex- 
tinction corrections used to calculate the Paa-luminosity Lpaa. 
The arrows visualize the change for extinction corrections done 
at progressively longer wavelengths. ULIRGs in which Hel 
2.058 //m is detected are marked by a square. 



4.3. [Fell] Emission in ULIRGs: Characteristic of tlie 
UNER-phenomenon 7 

lRiekeetalJ (ll98nl) detected for the first time a near-infrared 
line of [Fell] at 1.644 yum ( - a'D-112) in an extra- 

galactic object, the starburst galaxy M 82. A few y ears later, 
by inv estigating the supernova remnant IC 443, Gra ham et alJ 
('198'/) proposed that this line could be a good indicator for 
shock activi t y in a stronomical objects. Mouri et al . ( Js990) and 
iMouriet alJlll993h pointed out a correlation between [Fell] 1.644 
and [OI]63oo ascribed to the fact that [Fell] 1 644 and [OI]63oo 
arise mostly in partially ionized zones which are created either 
by relatively fast shocks (100 kms"') or by X-ray radiation 
from an AGN. They established a dignostic diagram, [Fell] 1 644 
and [OI]63oo being normalized by the hydrogen recombina- 
tion li nes Bry and Ha, respectively (Fig. la in iMouri et alJ 
('1990') and Fig. |3] in this paper). As the hydrogen recombi- 
nation lines arise from fully ionized zones, the [FeII]i.644/Br7' 
and [OI]63oo/Ha ratios are a measure of the ratio of the vol- 
ume of the partially relative to the fully ionized zones. Figure|3] 
shows the two line ratios: such a diagram sketches nicely the 
increasing importance of shocks, or of X-ray excitation when 
going from star-formation to AGN activity. Several studies 
(e.g. Moorwood & Oliva 1988: G reenhouse et al. 1991, 1997; 
iMouri et aljll993nForhes & WardH'l 993HVan^ Riekel[i997l) 
argued that in the case of starburst galaxies shock activity 
caused by supernova remnants and superwinds is responsible 
for the observed [Fell] emission and [OI] emission. 

Our SOFI spectra cover most of the H and K band, in- 
cluding [Fell] at 1.644 and/or 1.257 /im for several targets. 
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Table 5. Mean flux ratios [OIlesoo/Ho- and [FeII]i 644/Bry for 
subsamples of ULIRGs and emission-line galaxies 



Object class 


log{[OI]630o/Ha} 


log{[FeII]i.644/Br7) 


(1) 


(2) 


(3) 


LINER-ULIRG 


-0.93±0.07 


0.54±0.15 


HII-ULIRG 


-1.16+0.17 


0.19±0.08 


Seyfert2-ULIRG 


-1.00±0.07 


0.40±0.17 


Hll-galaxy 


-1.85±0.15 


0.10±0.10 


Seyfert-galaxy 


-0.98±0.18 


0.47±0.21 


Composite object* 


-1.57+0.09 


0.19+0.12 



: Composite objects exhibit both starburst and seyfert activity. 



This motivated a study of the [Fell] properti es of ULIRGs 
in general for which we added ULIRGs from IVeilleux et all 
lll997bl Il999b>) in order to improve the statistics. Our sam- 
ple includes ULIRG s with both rOIIfisnn and Ha measurements 
( Veilleux et alJl99ifl997illl999bl:lDuc et alJl997tlKim et alJ 
In total our sample comprises 21 ULIRGs. We catego- 
rize the ULIRGs according to their optically classified spec- 
tral type as Starburst-, AGN- or LINER-UL IRGs follow in 



the cla ssification scheme established by Ve illeux & Osterbroc 
lll987l) . In addition, we include in Fig.|3]for comparison, some 
lower luminosity Starbursts, Seyfert and composite galaxies. 
Representatives of pure shock- and photoionization excitation, 
the supernova remnant IC 443 and RCW 103 as well the Orion 
nebulae are added ^ . 

Since integrated ULIRG spectra represent a heterogeneous 
mixture of different astronomical objects, (e.g., Hll-regions, 
SNRs, PDRs, winds and AGNs) it is not surprising to find them 
between the loci of pure Hll-regions and pure shock excita- 
tion. Starburst-ULIRGs and LINER-ULIRGs populate diff^er- 
ent regions of the diagram. On average (see Table [S] and in 
Fig. |3] the enhanced symbols with error bars), both flux ra- 
tios, [FeII]i.644/Br7' and [OI]63oo/Ha, are higher in LINER- 
ULIRGs than in starburst-ULIRGs. The latter also have some- 
what higher ratios when compared to lower luminosity star- 
bursts. AGN-ULIRGs populate a region of high ratios, similar 
to lower luminosity Seyferts. The complex origin of the iron 
emission makes clear statements difficult, in particular for in- 
dividual sources. 

In what follows we describe a scenario which could be 
a plausible explanation of our observations: All starburst- 
ULIRGs and even more strongly the LINER-ULIRGs exhibit 
an elevated level of supernova- or superwind-related shock 
activity when compared to normal (lower luminosity) star- 
bursts. A number of effects may contribute to this behaviour: 
lAlonso-Herrero et alJ ( I2OO0I) . among others, suggested that the 



^ D ata of reference obj ects a re taken from: iMoQIffiflp d & Oliv^ 
jl988h: lArmus et alj ^1989^: IVeron-Cettv & Weroi^ 0986); 
KawaraetalJ h988h : Uitken et alJ <198lh : IPrenclJ <198Cl): 



[Johansson (198^ "Rieke et al.' ('1980); "Peimbert c& Spinrad 
(1970); Rieke & Lebpfsky. (198 1 ): Boksenberg et a l. (1975); Welch 



8^ 



Graham et alj 



197C); 
1987); 



Wilson et aLni97^lLester et al]"!^^ [Heck man et al 
CohenI <1983h; L owe et all jl979») nBalick et ali ti980) 



Jl987^: 



Essen & Kirshner ( 198( 
1983) 



lOliva et al 
( 19971 , 



Leibowitz & Danziaer ( 1983); Veilleux et al. 
Veilleux et aLll999a) ; .Duc et aLl 1997.) : . Kim et al..tl99^ ) 



198a 

1999bt) 



elevated shock activity (with respect to Hll-region activity) 
could be due to an aging starburst. The hot ionizing stars live 
only for about 10 million years while supernova and wind ac- 
tivity sets in after a few million years and persists longer The 
detection of the Hel 2.058 jim line in some starburst-ULIRGs 
but no LINER-ULIRGs is indeed consistent with this model. 
Furthermore, the measured equivalent width of the hydrogen 
recombination lines being an indicator of the age of the burst 
tends in LINER-ULIRGs to be smaller than in HII-ULIRGs 
Applying on our sample Fig. 7 of lAlonso-Herrero et al I teoool) 
where log{[FeII]i 644/Br')'}^ versus logfK/Bry) is plotted, we 
would estimate burst ages between 8-12 Myrs for the LINER- 
ULIRGs (assuming a burst duration of 1 Myr). Additionally, 
diff'erential extinction could also play a role: a larger scale 
wind will be less obscured than the embedded star forming re- 
gions. .Goldader et al.. (.1995.) invoked this scenario to explain 
the absence of a deficit in H2 emission in ULIRGs showing a 
Bry deficit. Qualitatively, such an effect may be indicated by 
sources having shock-like [OI]630o/Hq' ratios (from the outer 
wind) but more Hll-region like [Fell] 1.644 /Bry (from the more 
obscured starburst). Figure|3]may indeed contain such sources. 

AGN-ULIRGs populate a similar region in Fig. |3] as do 
lower luminosity Seyferts. This likely means a strong contribu- 
tion of AGN X-rays and AGN winds to their [Fell] emission, 
but the considerable overlap with other categories also allows 
significant non-AGN contributions. It is unlikely that AGN 
dominate the [Fell] emission of the starburst- and LINER- 
ULIRGs. It would be difficult to consistently avoid detection 
of [SiVI] from the AGN in that case. Most of the observed 
[FeII]/[SiVI] limits of the starburst-ULIRGs are higher than the 
ratios for bona fide Seyferts. This indicates that [Fell] emission 
in starbursts-ULIRGs and LINER-ULIRGs likely originates in 
a starburst superwind rather than in an AGN. 

4.4. Notes on individual galaxies 

4.4.1. F001 83-7111: Extreme [Fell] Emission 

F00183-71 1 1 is the most distant (z = 0.32) and most luminous 
object ( Li R = 10'^'^^ L(T|) in the ULIRG-sample observed with 



SOFT Ar mus et"aD ( 1 19891) classi fied this galaxy i n the optical 
as a LINER. In the mid-infrared. iTran et alJ (I2OOI) found that 
this object displays extraordinary absorption features of sili- 
cates and ices and interpreted it as a hea vily obscured AGN . 
Recently, Spitzer IRS observations by ISpoon et alJ ( l2004h 
strengthened this interpretation. Unfortunately, the prime AGN 
indicator [SiVI] is redshifted outside of the K band and cannot 
be used. We do not see a broad Pa/3 component. In this ULIRG, 
we measured the greatest intrinsic line width of [Fell] lines in 
our sample, ~ 900 km s" ' and ~ 1 000 km s" ' , for [Fell] 1 257 and 
[Fell] 1.644, respectively. The flux ratio of [FeII]i.644/Bry - 0.94 
(from measured Pa/3) is higher than in normal ULIRGs and 
similar to the 'near- ULIRG' NGC 6240. Both an AGN and a 
shock origin of this emission are possible, the large intrinsic 
half width of the iron lines, larger than Pa/?, may favour shock 



^ Following Simpson et al. ( 1996), this flux ratio could be used as a 
rough estimator for the age of the starburst. 
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Fig. 3. Logarithm of the flux ratio [Fell] 1644/617 versus the 
logarithm of the flux ratio [OIlesoo/Hor both for ULIRGs and 
lower luminous emission-line galaxies. The enhanced symbols 
with error bars represents the average of each ULIRG sub- 
sample. The dotted line shows a schematic mix of the Orion- 
nebula (Hll-region) and the supernova remnants IC 443 and 
RCW 103 (shocked sources). 

activity. The observed luminosity of both lines combined is 
3.7x10** Lq. For comparision, in NGC 6240 the observed lumi- 
nosity of these two Unes combined is 4.9 x 10^ Lq. The extinc- 
tion to the [Fell] emitting region can be estimated using the two 
iron lines arising from the same upper level, with an intrinsic 
flux r atio [Fell] 1.257/ [Fell] 1.644 of 1.36 jNussbaumer & StorevI 
1988'). We obtain A// = 1.45 mag using the extinction curve of 
Rieke & Lebofskv ( 1985), and Ah =1.2 mag using Koorneef 
ll983h . Assuming a foreground absorber we estimate an ex- 
tinction corrected luminosity of both lines combined of 2 x 10^ 
L0 and 1.5 x 10^ L0 for the two extinction laws, respectively. 
The extre me [Fell] luminosity and the unus ual mid-infrared 
spectrum jTran et alJ200lUSpoon et alJ2004 motivate further 
studies of F00183-7111 as an atypical ULIRG. 

4.4.2. F23578-5307: The galaxy with one of the 
highest H2 v= 1-0 S(l) luminosities 

The H2 lines of the rotation-vibration transition 1-0 S(...) domi- 
nate the spectrum of F23578-5307 (z = 0.125; LiR = 1.3xl0'2 
Lq) and, atypical for ULIRGs, Paff is not the strongest line 
in the near-infrared spectrum. Very few ULIRGs have simi- 
lar characteristics, the best known example being the 'near- 
ULIRG' NGC6240 which has long been f amous as the ^galax y 
with the highest luminosity in H2 emission (Joseph et al. 1984). 
Van de r Werf et al. (1993), Tecza et aL ( 2000) and Lutz et a l. 
( l2003b suggest that large scale shocks are responsible for the 
strong H2 emission in NGC 6240. Their line imaging in- 
dicates that the H2 emission in NGC 6240 does not arise 
from the two nuclei proper but from an extended region cen- 
tered between these two nuclei. Van der Werf et alJ ( Il993l) de- 
termined an observed luminosity of ~ 8 x 10^ Lq for the 



H2 v= 1 -0 S(l) line. In F23578-5307, the observed luminos- 
ity of the H2 v=l-OS(l) line is 3.1 x 10^ Lq, a factor 3.5 
more than for NGC 6240, and weU above flie ULIRG IRAS 
F09039-H0503 (1.1 x 10** Lrr.. IVeiUeux etail (Il997bh l. wifli 
one of the strongest H2 emission lines relative to Paa in the 
ULIRG population. The equivalent width of 42. 7A observed 
for H2 V = 1 -0 S(l) in F23578 -5307 is we ll above the - 20 A 
typical for ULIRGs ( Goldader etai]l995|) but short of the 7 1 A 
of NGC 6240. This confirms that the large H2 luminosity re- 
flects an unusual type of activity and is not just a consequence 
of scaling up a normal ULIRG spectrum. 

Further studies of these systems may help determining 
the origin of such extreme H2 luminosities. Is the situation 
for NGC6240, i.e. strong emission between two nuclei of a 
merging syste m, typical? Im ages of F23578-5307 taken at 
NTT/SOFI (see lRigODOuiouet aL.1999,) ) and the VLT/ISAAC 
(Tacconi, priv. communication) under good (0.8") and ex- 
cellent seeing (0.35") conditions, respectively, show that this 
ULIRG has a single nucleus. So far, there are no hints that this 
extraordinary ULIRG has a similar morphology as NGC 6240 
where the projected seperation of the two nuclei is 0.75 kpcs. 
At the distance of IRAS F23578-5307 this would correspond 
to X 0.4", allowing to barely resolve the two nuclei. 

5. Conclusions 

We have obtained moderate resolution near-infrared spec- 
troscopy of a sample of ULIRGs previously studied in the mid- 
infrared with ISO. A first goal was the search for near-infrared 
indications of AGN activity, namely broad components to the 
hydrogen recombination lines Paa or Pa/3, and detection of the 
coronal line [SiVI] 1 .962 jum. Most spectra do not show these 
indicators and appear starburst-like. In none of the 24 ULIRGs 
do we detect a broad velocity component of the hydrogen re- 
combination lines. In IRAS 041 14E we have evidence for a 
[SiVI] line with a luminosity of 1 .4 x 1 0* Lq (no extinction cor- 
rection applied). The agreement of near- and mid-infrared clas- 
sification is good, in line with previous suggestions that AGN 
in ULIRGs are often visible over a wide wavelength range, and 
with the majority of ULIRGs being predominantly starburst- 
powered. 

We have re visited the 'Bry-deficit' in ULIRGs found by 
iGoldader et"ai] (il995. . 1997b.) . Using the Une pair Pao', Ha for 
extinction correction, we confirm their result by finding a Paa- 
deficit. However, when using the two near-infrared lines Paa 
and Bry, the inferred extinction is significantly higher and no 
deficit of ionizing luminosity is found on average. This is con- 
sistent with most of the luminosity in these sources being pow- 
ered by star formation at Ay ~ 10-20 mag. 

Many ULIRGs are found to be strong emitters of [Fell]. 
ULIRGs that are optically classified as LINER or starburst also 
populate different regions in a diagnostic diagram using the ra- 
tios [FeII]i 644/Br7 and [OI]63oo/Ha. These observations are 
consistent with LINERs in infrared galaxies being due to super- 
nova and superwind related shocks. Age efifects and/or different 
location and obscuration of the wind and starburst regions may 
contribute to making a particular source a starburst or LINER 
type. 
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Fig. A.l. Kj-band image of IRAS 19458 plotted in grey scale 
and contour forms with a size of ~ 13" x 14". In an aperture of 
4 arcseconds, the magnitude in K, is 13.1 mag (Vega). Angular 
scale on the sky is indicated by the bar which represents 5 kpc. 
North is up, east to the left. 



Two unusual objects were identified in the sample. 
IRAS F00183-7111 which ISO observation suggest to be an 
obscured AGN exhibits extremely wide and strong lines of 
[Fell]. IRAS F23578-5307 is one of the most luminous in- 
frared galaxies in H2 v= 1 -0 S(l) according to our knowledge. 
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Appendix A: K band Image of IRAS1 9458+0944 

The 5 minute Kj-band image taken at NTT/SOFI during the 
second run tFig. lA.ll i shows a single system which seems to be 
almost relaxed. In an aperture of 4 arcseconds, the magnitude 
in K, is 13.1 mag (Vega). 
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